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Abstract 

Using quantum maps we study the accuracy of semiclassical trace formulas. The 
role of chaos in improving the semiclassical accuracy, in some systems, is demonstrated 
quantitatively. However, our study of the standard map cautions that this may not be 
most general. While studying a sawtooth map we demonstrate the rather remarkable 
fact that at the level of the time one trace even in the presence of fixed points on singu- 
larities the trace formula may be exact, and in any case has no logarithmic divergences 
observed for the quantum bakers map. As a byproduct we introduce fantastic periodic 
curves akin to curlicues. 
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1 Introduction 



The trace formula of Gutzwiller [1] provides a window into the semiclassical treatment 
of chaotic systems and is hence of considerable importance in many areas of physics and 
chemistry. Yet the understanding of this formula has been hampered by the complexity of 
the classical task involved, and often the semiclassical effort is greater than that for a full 
quantum solution which makes the semiclassical methods either only of academic interest or 
as tools for explaining features of a full quantum solution. In this situation simple models 
have played a key role just as has been the case with the understanding of classical chaos, 
as evidenced for example by the study of the logistic map and the smale horseshoe. 

Although the trace formula (and its relatives) are important, the accuracy of these have 
not been rigorously established. The trace formula being the first term in an asymptotic 
expansion, has many approximations with apparently uncontrolled errors. The very many 
applications of the stationary phase approximation, on functional and on ordinary integrals 
obscures the nature of the errors. One ongoing issue is the time at which we can expect 
the various semiclassical approximations to go bad with a candidate the so called log time 
which is of the order of — log(h) / A, applicable for chaotic systems with a Lyapunov exponent 
A, being called too pessimistic [2]. We note that although the semiclassical formula for the 
time evolution (due to Van Vleck) has been essentially known almost since the beginning of 
quantum mechanics the nature of the errors is still a matter of interest, especially since the 
trace formula depends on a Laplace transform of the Van Vleck approximation [1]. 

The quantum bakers map has been studied as a model of quantum chaos as indeed it 
has many generic features of quantized chaotic systems [3,4,5], however a rigorous study of 
some of the simplest traces of the quantum problem [5] showed logarithmic divergence from 
the expected Gutzwiller trace formula. The genesis of this disturbing fact was attributed 
to a fixed point on the line of discontinuity (corners) of the classical map. A baker without 
corners was studied [6], and it was found that there were no logarithmic divergences anymore. 
However we will in this note display a map, the sawtooth map, which has a fixed point on 
the line of discontinuity, but which has no logarithmic divergences. In fact we will show that 
in certain circumstances the semiclassical trace can be exactly the quantum trace. 

In pursuance of certain results in [5], we use the existing models of quantum chaos to 
study the accuracy of the simplest trace of the quantum problem. To this end we will in 
particular use the much studied quantum standard map and the much less studied quantum 
sawtooth map both in the settings of the phase space as a cylinder and as a torus. These 
being maps derived from time periodic systems we study simply the trace of the relevant 
unitary Floquet operators. In section 2 the cylindrical phase space is used, while in section 
3 the toral phase space is used, we end with some comments in section 4. In the Appendix 
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we provide a self contained derivation of the "trace formula" for a class of maps on the torus 
that includes the maps studied below and all the so called kicked systems. While this is in 
principle known, it may be useful to recollect it here, as the other published proof of Tabor 
[17] is for the case of maps on the entire plane, and is hence quite different. Those with little 
exposure to the trace formula may peruse the Appendix first. 

2 The Cylindrical Phase Space 

Consider the well known standard map (q n +i = q n + Tp n ,p n+ i = p n + K sm(q n+ i)) defined 
on the cylinder [0, 2n) x (— oo, oo) with periodic boundary conditions in q. The first iterative 
equation is them to be evaluated with a modulo 2n rule. A simple change of variables from 
p to Tp implies that the relevant parameter is just one in number and is k = KT. This 
area preserving map and its quantization have been studied in extensive detail as a model of 
chaos and its quantum manifestations and led to such important discoveries as the quantum 
suppression of chaotic diffusion. Some of these are reviewed in [7]. However its use in 
understanding the Gutzwiller trace formula has been limited. The work of Tabor [17] where 
he derived the relevant trace formula for quantum maps on the plane did use the standard 
map but no rigorous or extensive study has yet been attempted. The quantization of the 
standard map is described by an infinite dimensional unitary matrix with matrix elements 
in the momentum representation being 

U nm = {-i)" J n - m {k / Tft) exp (-im 2 aT/2) . (1) 

Here J v (z) is a Bessel function. As has been noted earlier there are two relevant parameters 
quantum mechanically and we can take T = 1 while we do not have the freedom to choose 
H — 1 anymore. The trace of the quantum map is then given by 

oo 

Tr(U) = J (k/h) Y z~ lhn2 ' 2 = Jo(k/h)e 3 (0\-h/2n) (2) 

n=— oo 

where 6 3 is a Jacobi theta function. An application of the Poisson formula enables us to 
rewrite the above as 

Phi °° 

Tr{U) = \ -rJo{k/h) Y exp(27r 2 m 2 /^) (3) 

Strictly speaking this sum need not exist and may without further regularization be infinite. 
This is not due to any peculiarities of maps, and in fact this is the rule in quantum mechanics, 
the trace of the propagator could be infinite in real time. However this will not worry us 
further as the semiclassical trace will also diverge exactly, allowing us to make meaningful 
relative error estimates. 
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Now we turn to a semiclassical evaluation of the quantum trace. We wish to then write 
the sum in terms of the period one or fixed points of the classical map. The Lagrangian 
generating the standard map on the cylinder is 

L(q n ,q n+1 ,l n ) = (q n+1 -q n + 2ixl n ) 2 /2 -kcos(q n+1 ) (4) 

as dL/dq n = —p n and dL/dq n+ i = p n +i- l n is the winding number in the q direction or 2nl n 
is subtracted by the modulo operation performed at time n to keep the map on the cylinder. 
The fixed points are determined by the condition (q n+ i = q n ,p n +i = Pn), and are therefore 
simply given by (g = 0,p = 2nl) and (q — ir,p — 2nl), where I runs over the integers. We 
note that the number of fixed points are infinite and that their positions do not depend on 
the classical chaos parameter k. The denominator in the Gutzwiller trace formula (or its 
versions for the case of quantum maps, they are essentially same) has the factor ^Tr( J) — 2, 
where J is the classical stability matrix associated with the periodic orbit. In the case of 
standard map fixed points this is either \f~k~ or y/—k depending on whether q = or q = tt. 
The semiclassical approximation is of the form 



Tr(U) sc = £ exp(iS/h)/ y (Tr(J) — 2) (5) 
fixed points 

where S is the classical action and since we are dealing with discrete time systems with the 
time interval being taken as unity this is identical to the Lagrangian written above. We get 
for the standard map on the cylinder the following semiclassical approximation for the trace 
(see Appendix): 

[~A k IT °° 

Tr(U) sc = J^ cos(| - J) £ exp(27r 2 *Z 2 M) (6) 

' l=—oo 

This then has to be compared with the Eq. (3). The asymptotic form of Jo(x) for large 
arguments is [8] 

/ 2 7T 7T 

J (x) ~ \ — [cos(x - -) P(x) - sin(x - -) Q(x)\ 
V nx 4 4 

where 

P{x) = 1 - 9/(128a; 2 ) + . . . and Q(x) = -l/(Bx) + ... 

where the neglected terms are at most of order x~ 3 . We then derive that 

Tl (U) ^.k, ,k ix x ^,k^ 

P (7-)-tMT--r)Q(T)- (7) 



Even this simple result is instructive. It indicates spectacular failures of the trace formula 
when the ratio k/h passes through 7r/4 + (2m + l)ir/2, where m is a positive integer. At 
these points the semiclassical trace becomes zero and the exact quantum trace is also small 
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indicating that such a failure will not be serious for the evaluation of the spectra. But it 
does suggest that the ratio of the quantum to the semiclassical need not be bounded in the 
entire range of the parameters. This divergence of the ratio is thus not due to any coalescing 
stationary points or bifurcations in general; as we have already noted the fixed points are 
independent of the chaos parameter k. Second, it is clear that for a larger k or chaos the 
semiclassical trace is going to be a better approximation. The authors in [2] have observed 
this, and we see that the simplest quantum trace provides a quantitative evidence of this. To 
compare the absolute error rather than the relative, we have to deal with the Jacobi theta 
function which occurs identically both in the quantum and in the semiclassical traces, and 
as this can depend critically on the rationality of n/% we do not pursue this here, except 
note that in unbounded phase spaces such as the plane and the cylinder the absolute value 
of the traces might actually diverge. 

The case of the standard map has been extensively studied but not so the quantum 
mechanics of the sawtooth map, which is in many ways simpler because it is a piecewise 
linear map for which it is easy to extract classical information. Consider the map (? n +i — 
q n + Tp n ,p n+ i = p n + K(q n+ i — tt)) defined on the cylinder [0, 2tt) x (—00,00), which is 
completely chaotic for positive K. This can be considered as a simple harmonically kicked 
particle in one dimension. We may view the system as an otherwise free particle on a unit 
circle being subjected to a kick whose strength is proportional to the physical angle and which 
has a single discontinuity when the angle is zero. The trace of the quantum propagator may 
be written down just as in the case of the standard map (T = l)with the result that 



where Erf is the error function. 

The semiclassical trace is also found quite easily, because the fixed points are at (q = 
n, p = 2irl), I being any integer. The classical action is similar to that of the standard map 
with the potential term being replaced by a periodicised quadratic function. Thus we get 
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sc 




(9) 



and the ratio of the quantum to semiclassical simply becomes 




(10) 



Inserting the asymptotic form for the error function we get the relative error to be 




(11) 
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The neglected terms are at most of order (h/K) 3 / 2 . 

In this map also we see the seemingly general result that the larger the chaos (larger K) 
the semiclassical accuracy is better. The precise way in which this occurs for this map is 
like K~ l l 2 . The difference between this map and the standard map is that the semiclassical 
approximation seems more well behaved and the sudden decline in the accuracy in the 
case of the standard map, due to the mismatch between the zeros of the quantum and the 
semiclassical is absent here. However in the regions where the semiclassical (and quantum) 
traces are not very close to zero, the standard map would perform better as the errors are 
of order h/k while for the sawtooth map the errors are of order (h/k) 1 / 2 . We will see that 
surprisingly many of our observations will remain valid in the toral setting as well. 

3 The Toral Phase Space 

Two dimensional maps on the torus (periodic boundary conditions in both position and 
momentum) have been studied extensively, from the standard map to the cat map of Arnold. 
Their quantizations have also been investigated by many [3,7,9], as this is a legitimate model 
for a bound Hamiltonian system with two degrees of freedom, the minimum needed for chaos. 
Unlike the case of the cylindrical phase space where scattering may take place, sheer topology 
now prohibits this on the torus. The act of taking periodic boundary conditions on the 
momenta is reminiscent of a similar procedure in condensed matter physics which has also 
provided an interesting two dimensional map, namely the kicked Harper model. Quantum 
mechanics now is on a finite dimensional Hilbert space and the we take the torus to be a unit 
torus (the dissected square is of unit sides). An alternate view of toral quantum maps is as 
the case of "resonance" of the cylindrical map in which case the quantum mechanics can be 
exactly reduced to that of a torus. Position and momentum take discrete values which we 
take to be n/N, where N is the dimension of the Hilbert space and n = . . . , N — 1. Here 
N is required to be an even integer and is the inverse of the Planck constant. Therefore the 
classical limit is iV — > oo. 

In this section we will first deal with the sawtooth map on the torus. We will write 
the map a little differently for convenience as (q n+ i = (q n + p„)modl, q n+ i = (p n + 
i^g n+ i)mod 1.). The quantum map or propagator is given by 

U nm = exp(i7rKm 2 /N) exp(ivr(n - mf/N). (12) 

v N 

This then is the quantum equivalent of the classical map restricted on a torus. It is a finite 
unitary matrix as all quantum maps on compact phase spaces are. Clearly it is one of the 
simplest quantum objects, having elementary functions as matrix elements. Yet we find that 
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it is rich enough to reflect many essential features of the effect of classical chaos on quantum 
systems. For more properties of a very similar quantum map please refer [10]. 
The classical Lagrangian (or action) for this map is 

L(q n ,q n+1 ,l n ,m n ) = (q n+1 -q n + /„) 2 /2 + K(q n+1 (modl)) 2 /2 - m n q n+1 , (13) 

where the integers l n and m n are winding numbers in the q and p directions respectively. 
The fixed points are given by (q = m/K,p = 0) where m — 0, 1, 2, . . . , [K\. The fixed point 
at the origin is in fact on the line of discontinuity, and is quite like the bakers map and its 
quantization. In that case it was noted that such singular periodic orbits, "half hyperbolic 
" , could give rise to unusual semiclassical behaviour and that the quantum and semiclassical 
could diverge semiclassically as \og(h~ 1 ). A bakers map was constructed to explicitly avoid 
this situation and a recovery of standard expectations was realized [6]. However if we take the 
above quantization as valid (since it follows canonical quantization closely [11]) it will turn 
out that the fixed point on the discontinuity is benign in the case of the sawtooth map, in 
that it does not lead to semiclassically divergent terms. This may suggest that the unusual 
nature of the bakers map originates either from the unusualness of the very quantization 
procedure adopted or is peculiar to the bakers map. 

Let us look at this case more closely, and take < K < 1, so that there is a single fixed 
point at the origin which is also a point of discontinuity. A direct application of the Gutzwiller 
formula will be incorrect in this case, but the rectification is obvious. The stationary phase 
integral must not extend over the entire real line, but must be restricted to [0, oo). Thus the 
semiclassical approximation is then the approximation (quantum ~ semiclassical) 



An estimate due to Wilton which we give below bounds the difference as N — > oo, but 
this is not a good estimate of the errors (for < K < 1), as these are of the order of 
the sums themselves. An application of the Euler-Maclaurin summation rule will give us 
the above approximation quite simply. But more important is the observation that it is 
legitimate to replace summations with integrals even in the presence of fixed points sitting 
on discontinuities, a rule that was flagrantly violated in the case of the bakers map. 

Let the errors be e(K, N) = \Tr(U) — Tr([/) sc |, that is the absolute value of the difference 
between the right hand and the left hand side of the Eq. (14). We now observe some 
remarkable number theoretic properties of these errors. First we note the surprising fact that 
if N is a multiple of 4, e(l/2, N) = 0. In this case not only does the semiclassical and quantum 
converge in the classical limit, they are exactly the same (at the level of the trace), despite 
the presence of the fixed point on the discontinuity! We have noted before that the main 




g-jyr/4 N-l 



1 



(14) 
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feature which makes the sawtooth map generic, in the sense that the semiclassical is only 
an approximation, and in the sense of eigenfunctions scarring and spectral random matrix 
properties is that the incomplete Gauss sums that arise while compounding propagators do 
not reduce to simple forms [10], unlike in the case of the cat maps (K an integer) where we get 
complete Gauss sums and has been investigated by Hannay and Berry [9] and the extensive 
number theoretic properties were explored by Keating [12]. Thus it is a bit interesting that 
an incomplete Gauss sum can be done exactly and is the semiclassical approximation. 
Consider the sum 

N-l 

a N =Y, expOn 2 /2iV), (15) 



n=0 

,2 



which is incomplete because the period of s n = exp(i7rn /2N) is 2N and the sum is only 
over the first N terms. First we note that Sjy/2+ m = s N / 2 - m exp(inm), thus terms with odd 
n do not matter as they cancel out of the sum exactly. Let 

2N-1 
n=0 

which is a complete Gauss sum [9]. But also note that if N is a multiple of 4, 

s n+N = (-l) n s n . 



This then implies that 



therefore 



Tr(U) = e- tn/4 a N /VN = l/y/2. (16) 



If N is not divisible by 4, then e(l/2, iV) is no longer zero and the semiclassical result is an 
approximation. Thus there are two distinct classes when K = 1/2 as far as the behaviour of 
the error is concerned, This can be generalized and we observe that if K = p/q, is a rational 
number, with relatively prime integers p and q, then there are as many classes as there are 
possible remainders of even N, on division by q. Within each class the error monotonically 
goes to zero in the classical limit, N oo. We have observed numerically that for large N 
to an excellent approximation 

e(K,N) = C{K)N- 1/2 . (17) 

Fig. 1(a) shows the case K = 2/5, when there are five classes. The constant C(K) 
depends on both K and the class, but we have observed that C{K) — C(l — K). If is 
irrational, the number of classes is infinite, there is no monotonic decrease of error, but the 
general trend is to decrease, as is evident from Fig. 1(b), where we have taken K to be the 
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Figure 1: The error as a function of inverse Planck's constant N, for the sawtooth map 
with (a, left) K = 2/5; (b, right) K = 7. 

most irrational number 7 = (\/5 — 1)/2, the golden mean. The number theoretical properties 
of the error is of course a reflection of the sum itself, and a fuller study is thus desirable. 

We however shall now turn to the case when K is larger than one, and there are a 
number of fixed points. If K is much larger than unity we need not correct for the corner 
fixed point as the errors incurred will be of the order K^ 1 / 2 , which is anyway the order of the 
total error according to Hardy, Littlewood and Wilton, as quoted in [13]. The semiclassical 
approximation is now described by 



Here [K] is the integer value of K. Note that this approximation clearly demonstrates 
the conjugacy between the Planck constant (1/iV) and the chaos parameter K. It is a good 
approximation if K is fixed and N tends to infinity, namely the classical limit, but also in the 
opposite limit of extreme "chaos", and small quantum system. For example when K tends 
to infinity for N = 2, a two level system. Undoubtedly, this approach to two level systems is 
to say the least inadvisable, as the quantum problem is much easier to solve. Nevertheless it 
throws up some interesting questions, because what would one mean by extreme chaos in a 
two level system? Quantum chaos for legitimacy in using the epithet chaos has always relied 
on the classical limit, and in the extreme quantum limit the concept becomes nebulous. 
The above approximation is derived from the Gutzwiller formula, using the [K] + 1 fixed 




e -in/4 N-l 



1 



J2ex-P(-inNn 2 /K). 



(18) 
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point actions and the stability matrix. The fixed point at the origin is included wrongly 
in the formula and we have discussed it separately in the case < K < 1 above. Now we 
use an estimate of Wilton, quoted in [13], to bound the error in this approximation. This 
estimate is the inequality 



\Tr(U) - Tr(U) t 



< 



N 



+ 



K 



(19) 



where a and b are constants of order one whose values are slightly different from the quoted 
ones, as the sawtooth sum is from to N — 1, rather than 1 to N, but this is a trivial 
complication. The essential conclusion is that the error is bounded in N, and corroborates 
our conclusion that there is no log(iV) problems in the sawtooth map despite the presence 
of the fixed point on the singularity. Also we once again notice that a larger K (larger 
chaos) implies a better semiclassical accuracy. In Fig. 2 we show the results of a numerical 
computation of the error, which shows that the error goes roughly as K~ 1 / 2 , with the slope 
depending weakly on iV and K and rather strongly on the fractional part of K. Notice that 
at integer K (cat maps) the approximation is exact. We thus note that like its relative on 
the cylinder the torus map error is also of the order of K~ 1 / 2 . 




Log(K) 



Figure 2: The error as a function of the chaos parameter K, for the sawtooth map with 
N = 4000. 



We had earlier studied the sawtooth map on the torus with slightly different boundary 
conditions [10], and came to essentially similar conclusions, however in that case there was 
no fixed point on the discontinuity. When the fixed points moved towards the singularity, 
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at the exact point of contact, the map turned into the continuous cat map when the trace 
formula becomes exact. We now turn to a smooth map on the two torus, the standard map. 
Much studied in other contexts, it is largely unstudied from the view of trace formulas, 
although it should form a natural object of study, as it represents generic mixed systems, 
and one recognizes the difficulties of such systems apropos trace formulas. 

We will consider the map once again on the unit torus. The map is now (q n+ i = (q n + 
p n )modl, Pn+i = (Pn + k sin(27rg n+ i))modl). The Lagrangian (action) is 

L(q n ,q n+1 ,l n ,m n ) = (q n+1 - q n + Q 2 /2 - k cos(2irq n+1 ) /2ir - m n q n+1 . (20) 

The fixed points are easy to find and are at (q m = ^ arcsin(m//c),p m = 0), where m is an 
integer such that \m\ < k. Once again the location and number of the fixed points are very 
much dependent on the the chaos parameter k unlike the standard map on the cylinder. 
Note that there are two fixed points associated with each m. 

The quantum propagator is again given by (N even throughout) 

e -i7r/4 

U nm = -^=exp(-ikNcos(2irn/N))exp(iiT(n-m) 2 /N). (21) 
v N 



and the trace is the finite sum 

e -j7r/4 N-l 

Tr(U) = —^J2exp(-ikNcos(27rn/N)). (22) 

The semiclassical approximation must approximate this sum using the classical fixed points, 
their actions and their stabilities. First note that the quantum trace apart from the factor 
g-jTr/4 j g rea j_ rjj^ sawtooth ma p traces were Gaussian sums, and partial sums of these give 
rise to the "curlicues", that have been studied in detail earlier and find applications from 
optics to quantum mechanics [13]. The partial sums associated with the quantum trace of 
the standard map ought to be interesting mathematical objects in their own right. 

Here we present a preview of the intriguing class of periodic "curves" that this gives rise 
to. Consider the sum 

i 

a N (k,l) = ^exp(-ikN cos(2irn/N)). (23) 

n=0 

The quantum standard map trace is essentially <7jv(fc, iV). We plot in Fig. 3 the real, x(l), 
and imaginary, y(l), parts of a^(k, I) as an argand diagram. Strictly there are only as many 
points as terms in the partial sum, and as such do not define a curve, but the connections 
showing the motion of the points in "time" I is suggestive. Only the sum up to / = N/2 
is relevant as beyond this a reflection and periodic translation is the whole curve. Hence 
unlike curlicues, these are periodic curves. A plethora of cornu spirals interweave to conjure 
the extravagantly ornate features. Several "eyes" of cornu like spirals are connected in a 
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remarkable manner. An essential feature of these curves is their extreme susceptibility to 
small changes in k, especially at large N values. We cannot rule out the possibility of number 
theoretic features here as well, and indeed we observe marked difference between the cases 
when k is an integer and when it is not. 

The periodic nature of the curves are easily understood, as y(l + N/2) = —y(l) and 
x(l + N/2) = x(l) + A, where A = J2 n =o,N/2-i cos(—ikNcos(2nn/N)). Hence the period of 
the curves are 2A. As iV is increased for fixed k the curves tend to smoothen out more but 
with the cornu spiral structures intact. Also note that the quadratic approximation of the 
cosine term in the exponential are the culicues studied earlier [13]. 




-150 -100 -50 -40 -20 

K=10 K=10.1 




-40 -30 -20 -10 20 40 60 

K=10.4 K=10.9 



Figure 3: A panel of the partial sums a W00 (K, I) (Eq. 23) for various K. 
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A detailed study of these curves is however not pursued in this paper, and we turn once 
more to the central theme, namely the accuracy of semiclassical approximations. Define 

6 N (l) = N(Vk 2 - I 2 + /arcsin(//A;)). (24) 

Then a careful application of the trace formula yields the semiclassical approximation 



p -iTv/4 

Tr(U) sc = 



2tt 



2 sin(iV£; + tt/4) 4 sm(8 N (l) + tt/4) 



(25) 



which is real except for the factor e - " 1 "/ 4 , just as the exact quantum trace in Eq. (22). Here 
the fixed point at the origin has been treated separately as corresponding to each \m\ there 
are four fixed points except when m — 0, which gives rise to two fixed points. The range of the 
inverse sine function is [0,7r/2]. Note that this approximation can be derived independent 
of the trace formula, by first converting the sum into an infinite sum of integrals via the 
Poisson summation formula, and of these integrals retaining only those which contribute to 
the stationary phase, and evaluating these integrals in the quadratic approximation. Indeed 
the trace formula in the context of maps on the torus is derivable by this procedure and is 
hence considerably simpler than the derivation for general Hamiltonian flows (see Appendix). 

First we note a drawback of the semiclassical approximation. Near parameter values 
where bifurcations are going to create new periodic points the semiclassical approximation 
becomes poor. Uniform approximations may be made to smooth over these parameter values. 
For instance in the case of the standard map, whenever k passes an integer value it creates 
four new fixed points via tangent bifurcations, and the approximation to the trace becomes 
poor. This is illustrated in Fig. 4(a). However notice that the approximation is indeed quite 
good even just away from these values. We will below retain the approximation made above 
and keep away from points of bifurcations to study various aspects of the error. 

In Fig. 4(b) we show the absolute value of the error e = \Ty(U) — Tr(U) sc \ for a fixed 
value of k — 100.5, and a range of N. It is clear that the errors involved are not monotonically 
decreasing with N, but that there are fluctuations in the mean. The mean however very 
nearly follows N^ 1 behaviour. This should be contrasted with the more cleaner but slower 
N~ l l 2 behaviour of the sawtooth map. In general the errors incurred for the smooth standard 
map is smaller than those of the sawtooth map. We see that this is a property that is carried 
over from the maps on the cylinder. It is unclear whether the standard map fairs better due 
to its smoothness or due to some other property. This said, we note that the trace formula 
approximation behaves better with increasing chaos in the case of the sawtooth map than 
the standard map. 

At a fixed value of N the behaviour of the trace for increasing k is quite complex. No 
simple rule like, higher chaos implies better semiclassical accuracy emerges. First we note 
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Figure 4: The error for the standard map on the torus as a function of (a, left) k, with 
N = 1000; (b, right) inverse Planck's constant N, with k = 100.5 

that the ratios or the relative error fair quite poorly, with large abrupt deviations. Fig. 
5(a) shows the relative errors |e|/|Tr([/) sc ], but the range in the y axis has been lessened 
to highlight the predominant features. The error could become as high as 40 in the range 
of k shown. This is however easy to understand in the context of the Eq. (7) for the map 
on the cylinder, as due to the mismatch of the zeros of the semiclassical and the quantum 
traces. We expect the absolute value of the errors to fair better. Indeed when N is small, 
one observes an increase of accuracy with k, and this is a power law for small N. Thus we 
show in Fig. 5(b) the logarithm of the errors with the logarithm of k, and we observe the 
rule that e = Ck -1 , with 7 ~ .25, for the extreme quantum case N = 2. Here the k go in 
steps of unity with their fractional part being uniformly .5. 

However as N is increased we observe that this is no longer true, and that the errors 
tend to oscillate with A; in a fairly complex manner. The overall tendency seems to be in 
fact a deterioration of the accuracy with k. Fig. 8 shows the errors as a function of k 
in a large range. The interesting aspect is that the oscillation period seems to increase 
geometrically. Therefore there is a possibility that the errors scale with k, but in a complex 
manner. However increasing k even higher seems to render the errors an almost monotonic 
growth, and we cannot comment on the many details until further analytic estimates are 
made and computer generated errors are under control. 

An important difference between the semiclassics of the sawtooth map and the standard 
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Figure 5: (a, left) The relative error for the standard map, as a function of the chaos 
parameter k, N = 1000; (b, right) the absolute error as a function of k in the extreme 
quantum limit, N = 2. 

map is that because of the piecewise linear nature of the sawtooth map there was no need 
to truncate Taylor expansions about the periodic orbits, while in the case of the standard 
map this truncation is done up to second order. Higher order corrections can be worked out 
and have formed the subject of recent works [15,16], it should be interesting to compare the 
errors when such higher order terms are included. 

4 Summary 

Using maps as models of quantum chaos we have studied the simplest traces from the point of 
view of semiclassical approximations, hoping that this will reflect some more general aspects 
of errors involved in the Gutzwiller's periodic orbit sum or trace formula. We found that 
indeed the maps on the cylinder had controllable relative errors. We found that the smooth 
standard map faired better than the sawtooth map in the order of magnitude of the errors. 
We also noted that the mismatch of zeros between the semiclassical and the quantum can 
lead to deterioration of relative errors. In the case of toral maps we displayed a rather 
remarkable fact that even in the presence of fixed points on discontinuities the trace formula 
may be exact at the level of the time one trace. Otherwise we showed that more generally 
there is no log(fr) divergences in the sawtooth map, unlike the case of the quantum bakers 
map. We also noted that the error may actually decrease with increasing chaos, a rather 
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Figure 6: The error as function of A; in a wide range of k, for the standard map with N = 1000. 

counter intuitive fact, which has quantitative evidence in the case of the sawtooth map. 
However we caution that this cannot be a general principle as a study of the standard map 
has shown a much more complex behaviour of the error. It is indeed a matter of importance 
to study and contrast the two maps, and understand why the errors behave qualitatively 
differently. 

We have restricted our study to the time one trace and indeed it is natural and necessary 
that we study the higher traces. For instance to study the errors involved in the semiclassical 
evaluation of Tr(U 2 ) leads us into the problem of enumerating period two orbits, which can 
be fairly easily computed in the cases studied, but which are much larger in number than 
the fixed points. Preliminary results in this direction indicate that for the sawtooth map the 
error behaviour of the time two trace is similar to that of the trace, although the errors are 
larger. We have also displayed some remarkable curves akin to curlicues that arise out of a 
study of the trace of the standard map on the torus. 
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Appendix 



Consider the following map denned on the unit torus: 

(26) 



q j+1 = qj + pj - lj 



Pj+i = Pj - v'(qj+i)- m j 



j- 



where V(q) is a periodic "kicking" potential with unit periodicity, the prime indicates deriva- 
tive, and the integers rrij and lj are such that the map is restricted to the unit torus (winding 
numbers). They are the "modulo one" operations. As is well known the above map can be 
derived from a time dependent Hamiltonian system. Also there exists a generating function 
for the map on the torus given by 

S(qj, q j+ i] lj, mj) = ^(qj+i ~ % + lj) 2 ~ V(q j+1 ) - m jqj+1 (27) 

from which the map can be derived as pj = —dS/ dqj and pj + i = dS/ dqj + \ . 

The quantum map corresponding to the Eqs.(26) is given by the finite dimensional unitary 
matrix 

(n\U\n') = ^jt exp(^(n - n'f - 2niNV(^)). (28) 



This takes on the function of the propagator as it connects states after consecutive kicks. 
Here the representation is in the position basis and n and n' may take the integer values from 
to N — 1. N is itself the dimensionality of the Hilbert space and is also the inverse Planck 
constant, and is restricted to be even integers for preservation of toral boundary conditions. 
Semiclassics means the study of the unitary operator as N — > oo. 

The object of primary interest below is the trace of the powers of the propagator from 
which the spectrum may be derived by a Fourier transform. This can be written as 



-inT/4 N-l 1 2 • T , s 

7li — \ 7 — 1 



(29) 



with 



riT+i = ni. 



It is implied in the above that the outer sum is over the T variables rij, % — 1, 2, . . . , T. Using 
the Poisson summation formula we can rewrite the above as 



Tr(f/ T ) = e -™ T ' A N T l 2 ]T J € d Xl ...dx T 

ki=—oo ~ e 

exp j 27niV^ (kjXj + &j — x j x j+i — V( x j)) ) ? (30) 
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with 

xt+i = Xi, and kx+i = k±. 

The stationary phase approximation, assuming that N is large, gives the following T condi- 
tions 

x j+ i = 2xj - Xj-i - V'(xj) + kj\ j = 1, . . . ,T, x = x T , xi = x T +i, (31) 
with the further condition that 

< Xj < 1. 

These are precisely the equations that determine the period T orbits of the map. The integers 
kj are related to the integers lj and m 3 - of Eqs. (26) by the relation kj = — lj — mj-\. 
With the assumption that they uniquely determine a periodic orbit of period an integer 
fraction of T, we restrict the infinite sum to only those that relate to such a periodic orbit 
labeled below by (3. This is the first approximation and the second one is that we will extend 
the ranges of integration to the whole real line. These approximations together are exact 
in the case of the cat maps [9,11,12], and the third related approximation is that we Taylor 
expand the potential about the periodic orbits retaining up to the quadratic terms. This 
approximation is unnecessary for the piecewise linear sawtooth maps [10], but is necessary 
in general. 

Therefore we define new variables yj such that Xj = qj + yj and qj is the central periodic 
orbit. 

Tr(U T ) ~ e-^ 4 N T / 2 £ exp ( 2niN £ (- mj+1 + q) + q 3 k 3 - V{q 3 ))\ x 

P \ 3=1 J 

J^d yi ...dy T exp (27TiNj2(y^ - y m _ x - y*V" /2)j (32) 

The sum 

T 

(-QjQj+i + 9j + Qj k i - v (Qj)) ( 33 ) 

i=i 

can be identified with the action of the periodic orbit, as from Eq.(27) one gets 

Sp = E {<$ - QjQj+i ~ V( qj ) + qjkj + Zj/2), (34) 

3=1 

which is essentially the sum in Eq. (33). We can neglect the term /|/2, from the Lagrangian; 
this is legitimate as the action occurs as e 2mNS , and since iV is an even integer these terms 
do not matter. Note that the trace formula at T = 1 is easily obtained from the above and 
has been used in the main body of this paper. 
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The integral is evaluated by standard methods. 

/oo 
d T y expiiirNytApy) = 
-oo 



e i7rT/4 



(35) 



Here y denotes the T component vector {yj} and y l is its transpose. The real symmetric 
matrix Ap is given by 



(2- V"( qi ) -1 . . . . -1 \ 

-1 2 - V"(q 2 ) -10.. 



Ar = 



V 



-l 



o 



(36) 



-1 2 - V"(q T ) J 



up is a "Maslov like" index, and is the number of negative eigenvalues of the matrix Ap. 

The nontrivial problem of evaluating the determinant is already solved and the results 
are well known [14]. The determinant can be related to the stability of the periodic orbit. If 
the residue of the periodic orbit, which is simply related to the trace of the stability matrix, 
is Rp, then Rp = — ^Det(Ap). or if we denote the stability matrix eigenvalues by and 
A^, Det(Ap) = + -2. 

Thus finally we can write the asymptotic periodic orbit sum, or the trace formula, for 
the toral maps as 

^ ^ exp (2niNSp — iirvp/2) 



Tr(f/ J 



E 

P 



(37) 



The sum is over periodic orbits whose periods are integer fractions of T are labeled by /3. This 
is of the canonical Gutzwiller- Tabor form having in the exponent the action of the period T 
orbits, and the prefactor explicitly depending on the stability of these orbits. The piecewise 
linear maps such as the baker's map, or the sawtooth map are such that the Maslov like 
phases are zero and the only type of periodic orbits are those of the direct hyperbolic kind. 
Now however we can study the vast class of mixed systems, such as the well known standard 
map. It is then of interest to ask questions such as how long the semiclassical approximations 
are valid in such systems which one knows are generic. The quantum mechanics on the torus 
is exact and thus we expect that the derivation above will facilitate the investigation of the 
periodic orbit sum for mixed systems. In this note we have used the above formula only in 
the special case of T = 1 . 
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